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Abstract. We study the dynamic behaviour of concentrated colloidal hard 
spheres using Time Resolved Correlation, a light scattering technique that 
can detect the slow evolution of the dynamics in out-of-equilibrium systems. 
Surprisingly, equilibrium is reached a very long time after sample initialization, 
the non-stationary regime lasting up to three orders of magnitude more than the 
relaxation time of the system. Before reaching equilibrium, the system displays 
unusual aging behaviour. The intermediate scattering function decays faster than 
exponentially and its relaxation time evolves non-monotonically with sample age. 



1. Introduction 

In the last two decades, colloidal hard spheres have been widely studied as a model 
system that exhibits a transition to an arrested glassy state as particle volume fraction, 
<p, increases £Q. Above tp w 0.49, the dynamics dramatically slows down and exhibits 
a two-step relaxation, where the initial relaxation corresponds to the motion of the 
particles within the cage formed by their neighbours, while the final, a-relaxation 
is due to cage-escape processes El- Probably, the most popular experimental 
technique for investigating these dynamics has been dynamic light scattering, which 
has allowed the microscopic dynamics of samples at thermal equilibrium to be probed 
for volume fractions up to p ~ 0.55 3 . Light scattering experiments measure the 
time autocorrelation function of the fluctuations of the scattered intensity, g 2 (r) — 1, 
which is proportional to the squared intermediate scattering function 0]. Confocal 
microscopy was used in more recent works that focused on dynamic heterogeneity, 
identifying localized clusters of particles whose dynamics differ from the average |SJE]. 

By contrast, experiments probing the non-equilibrium or aging behaviour of 
concentrated samples are scarce [31 El E]- This is because conventional light 
scattering techniques require a considerable time averaging to get sufficient statistics. 
Clearly, extensive time averaging prevents aging dynamics to be precisely characterised 
in systems where the dynamics are not stationary. Microscopy experiments can in 
principle provide time-resolved information, but they have often very limited statistics, 
since they are forced to deal with rather large particles, for which the dynamics is too 
slow to be followed in the aging regime over several decades. 

In this work, we take advantage of recently-introduced dynamic light scattering 
techniques, multispeckle ^J^| and Time Resolved Correlation (TRC) to follow 
the slow evolution of the dynamics of concentrated suspensions of colloidal hard 
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Figure 1. Sketch of the experimental setup described in sectiori f2.ll 



spheres, from an initial aging regime up to times where thermal equilibrium is reached. 
We find that hard spheres display an aging behaviour that is quite different from 
the one observed in molecular supercooled liquids. Strikingly, thermal equilibrium 
is reached only for times that considerably exceeds the relation time of the system 
and an unusual 'compressed' exponential relaxation of (72 — 1 is observed in the aging 
regime. Additionally, we find that in the stationary regime the dynamics appears 
to be extremely sensitive to a small shear, an issue possibly relevant when devising 
experimental methods dealing with very concentrated samples. 



2. Experimental methods 

2.1. Experimental setup 

The experimental setup is shown in figure The light source is a frequency-doubled 
Nd:YAG laser whose vertically polarized beam impinges on a cylindric scattering cell, 
C. The cell temperature is set to T = 20 ± 0.05°C via water circulating in a copper 
cell holder, T. Light scattered within a volume of about 1 mm 3 at a scattering angle 
9 = 90° (which corresponds to wavevectors close to the first peak of the structure factor 
of our hard sphere system) is split by the beam splitter BS so as to illuminate a charge- 
coupled device (CCD) detector and to be collected by a single-mode, polarization- 
maintaining fiber optics. Two photomultipliers connected to a hardware correlator 
hosted in a personal computer (PC) allow the fiber optics signal to be analyzed to 
obtain g 2 ( T ) — 1 down to a lag time r = 50 nsec, for samples whose relaxation time 
does not exceed 100 sec and whose dynamics are stationary (here, 52 — 1 is obtained 
by extensive time-averaging). 

The CCD images of the speckle pattern generated by the sample are transferred 
to the PC via a frame grabber and stored on the hard drive for later processing via a 
software correlator. The quantity measured in our TRC experiments is the degree of 
correlation, c/(t, r), of the scattered intensity of two images separated by a lag r |12j : 

r(f ^- (I P (t)I p (t + r)) p 

(I P (t)) P (I P (t + r)) P W 

where I p {t) is the instantaneous value of the scattered intensity measured by the p-th 
pixel at time t, and the average is taken over all CCD pixels. Note that no time 
average is performed; thus, the degree of correlation depends on both the time lag 
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r between the two images analyzed and the time t at which the first image of the 
pair is taken. In systems with stationary dynamics a further average over t can be 
performed. The normalized auto-correlation function is then obtained 

92(r)-l = (c I (t,T)) u (2) 

similarly to conventional dynamic light scattering techniques. The minimum time lag 
for which c/ can be calculated, t ~ 20 ms, is limited by the speed of the CCD camera 
and the exposure time. By coupling DLS data obtained from the hardware correlator 
to TRC, we are able to measure simultaneously and on the same sample the intensity 
correlation function from r = 50 ns up to very large lag times, r ~ 10 4 s. 

2.2. Sample preparation 

The particles used in this work are poly-methylmethacrylate (PMMA) spheres of 
radius R = 140 nm, with a degree of polydispersity of about 10% that prevents 
crystallization in the experimental timescale. The particles are sterically stabilized by 
a thin layer of poly-12-hydroxystearic acid of approximately 10 nm thickness. They 
are suspended in a mixture of organic solvents (cis-decalin and tetralin) that closely 
matches their refractive index and reduces van der Waals interactions. 

Samples are prepared by centrifugation of a stock suspension for about 24 h at 
1000 g to obtain a dense amorphous sediment, assumed to be near random close 
packing, tp ~ 0.67, for a 10% polydispersity 13 . The clear supernatant is removed 
and a controlled mass of solvent is added to obtain the desired volume fraction, 
ip. Samples are then dispersed by vortexing and tumbling for about 6 h. Before 
starting measurements the sample is thermalized at 20°C for 15 minutes. In aging 
experiments, the waiting time, t w , is taken from the end of the dispersion, which 
corresponds to t w — 0. Since our sample is very polydisperse we cannot use the 
freezing volume fraction to precisely control volume fractions El, leading to 
rather large uncertainties on the calculated tp. Indeed, a comparison of our dynamic 
data with the literature suggests that we systematically underestimate volume 

fractions by 2 to 4%. This should be kept in mind when comparing the present results 
to the literature. 

3. Results 

The left panel of figure El shows the waiting time evolution of the degree of correlation 
ci(t W7 r) measured for a sample at ip — 0.545, for various time lags r ranging from 
2 to 2000 sec and corresponding to the a-relaxation of the system. At fixed lag, ci 
continuously increases with t w , demonstrating that changes of the sample configuration 
during a fixed time interval become smaller when i w increases. In other words, the 
dynamics slows down when the system gets older, a clear manifestation of aging |14| . 
The right panel shows a set of representative intensity correlation functions gi (t w , r)— 1 
obtained by averaging c/ over a short time window (a few hundreds of sec) , starting 
at various t w . To extract the characteristic time, r a (t w ), of the a-relaxation, we fit 
the curves by a stretched exponential: <?2 (^w , T ) — 1 = acxp [— (r/r Q ) p ]. The inset 
shows r Q vs. i w , indicating that the relaxation time increases linearly with sample 
age, a behaviour observed in many glassy systems [TJj and refereed to as 'simple 
aging'. Note however that r a /t w -C 1, thus ruling out the simple explanation stating 
that the relaxation time is set by the age, r Q r* t w . Surprisingly, we find that for all 
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Figure 2. Left: age dependence of the degree of correlation, cj(t w ,r), measured 
for a sample with ip = 0.545. From bottom to top, r = 2000, 1200, 1000, 800, 600, 
400, 300, 200, 120, 40, 2 sec. The steady growth of c 7 (t w ,r) with t w indicates 
aging. Right panel: time decay of ^(iwjT) — 1, obtained by averaging c/(f w ,r) 
over a time window of a few hundreds sec, for t w = 8400, 12400, 18400, and 
29400 sec (from left to right). The shape of the relaxation is well described by 
a 'compressed' exponential form, as shown by the full line which has a streching 
exponent p = 1.22. Inset: Linear dependence of the a-relaxation time on t w 
typical of a 'simple aging' behaviour. 



curves p s=s 1.3. This behaviour is in contrast with that of concentrated colloidal hard 
spheres at equilibrium, where p< 1 0. Intriguingly, it is similar to the 'compressed' 
exponential relaxation found in a variety of soft glassy materials out of equilibrium |15| . 

The sample shown in figure|21does not reach equilibrium on the experimental time 
scale. In order to investigate how a stationary state is eventually reached, we turn to 
a more diluted sample, ip = 0.505. Figure |21 shows the t w dependence of c/. A striking 
behaviour is observed. Initially, cj grows, as also observed in figure [21 However, 
after reaching a maximum, the degree of correlation then decreases with £ w , before 
eventually reaching a stationary regime. This implies that after the initial simple 
aging regime, the dynamics accelerates and finally becomes stationary. Moreover, for 
different r the maximum of ci occurs at different ages, implying that also the shape 
of (72 — 1 changes with time. 

The evolution of Q2 — 1 is shown in figure^ Initially the relaxation time increases 
and the correlation functions are compressed [p w 1.3). Then, r a goes back to 
smaller values. Concomitantly, the shape changes and becomes more stretched. For 
tw > 1.8 x 10 5 sec the dynamics becomes stationary. Thus, equilibrium is reached after 
more than 1000 relaxation times, a behaviour at odd with that of molecular glassy 
systems This suggests that motion on length scales much larger than the particle 
size might be involved in this process. At equilibrium, we measure the fast dynamics 
with the hardware correlator by averaging gi — 1 over 12 h. The inset of figure 0J shows 
the full (/2 — 1 obtained by pasting together hardware correlator and CCD data. A 
typical two-step decay with a final stretched relaxation (p w 0.6) is observed. 

Values of the t w dependence of r Q and p are shown in figurc[3for the two samples 
discussed above and for an additional intermediate volume fraction. For the lower tp 
(open circles), three regimes are observed: aging (r a oc t w and p ~ 1.3), acceleration 
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Figure 3. Degree of correlation as a function of t w for a sample with ip = 0.505. 
From bottom to top, r = 200, 140, 100, 70, 40, 30, 18, 10, 4, 2 sec. Before the 
x axis break, the scale is logarithmic, in order to better appreciate the initial 
evolution of the dynamics. The rectangles indicate the time windows over which 
cj is averaged to obtain g% — \ shown in figure IH 
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Figure 4. Intensity correlation functions calculated from the cj of figure El Each 
curve is labelled by the time window shown in figure 13 The initial aging regime 
(1 — > 2 — > 3) is followed by an acceleration of the dynamics (3 — » 4 — » 5) and 
finally by a stationary regime (5). Inset, CCD data (solid squares) and hardware 
correlator data (open squares) collected simultaneously in the stationary regime. 



of the dynamics (r Q and p decrease with t w ), stationary dynamics (r Q and p constant). 
As the volume fraction increases, the end of the aging regime shifts towards larger t w , 
eventually becoming too large to be detected. 

Our findings indicate that it takes a surprisingly long time for concentrated hard 
spheres to recover from a mechanical perturbation, such as the initial homogenization. 
This equilibration time may become prohibitively long when r a exceeds 1000 s. 
Moreover, ultra-slow equilibration is potentially relevant to several measuring schemes 
used to deal with glassy dynamics, where the sample cell is rotated to accumulate 
statistics, either in between short DLS runs when measuring short time dynamics 
('brute force' method JSj), or continuously to access the slow dynamics ('interleaved' 
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Figure 5. Age evolution of the relaxation time T a (left panel) and the stretching 
exponent p (right panel) for samples at ip = 0.545 (open triangles), tp = 0.54 (solid 
squares), and ip = 0.505 (open circles). The lines r a oc t w and p = 1.3 are guides 
to the eye. 
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Figure 6. Time evolution of the degree of correlation for a sample with 
<p = 0.496, after reaching equilibrium (from bottom to top, r = 200, 100, 30, 
2 sec). The sample is slightly rotated at three different times indicated by the 
arrows. Recovery after each rotation is very slow. 



[XT] or 'echo' JS] methods). Indeed we show the effect of a tiny rotation of the cell 
in figure |5| Initially, the degree of correlation is almost constant, since the sample is 
at equilibrium. When the sample is turned manually by the smallest possible amount 
(arrows in the figure), cj drops to 0, because the speckle pattern changes completely. 
Strikingly, it takes a very long time for c/ to grow back to its initial value. Moreover, 
the larger r the longer the recovery time. We note that at the smallest delay shown in 
figureEJ cj stabilizes in a few seconds, suggesting that the short time dynamics should 
not be significantly affected by the sample rotation. By contrast, on the longer time 
scales probed by the interleaved or echo methods, sample rotation may spuriously 
accelerate the dynamics. It would be interesting to also test these effects in setups 
where sample rotation is done by a stepper motor. 
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